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The dynamical properties of an extended Hubbard model, which is relevant to quarter-filled lay- 
ered organic molecular crystals, are analyzed. We have computed the dynamical charge correlation 
function, spectral density, and optical conductivity using Lanczos diagonalization and large- -/V tech- 
niques. As the ratio of the nearest-neighbour Coulomb repulsion, V , to the hopping integral, t, 
increases there is a transition from a metallic phase to a charge ordered phase. Dynamical prop- 
erties close to the ordering transition are found to differ from the ones expected in a conventional 
metal. Large-iV calculations display an enhancement of spectral weight at low frequencies as the 
system is driven closer to the charge ordering transition in agreement with Lanczos calculations. As 
V is increased the charge correlation function displays a plasmon-like mode which, for wavevectors 
close to (ti,h), increases in amplitude and softens as the charge ordering transition is approached. 
We propose that inelastic X-ray scattering be used to detect this mode. Large- N calculations predict 
superconductivity with d xy symmetry close to the ordering transition. We find that this is consistent 
with Lanczos diagonalisation calculations, on lattices of 20 sites, which find that the binding energy 
^ 1 of two holes becomes negative close to the charge ordering transition. 

PACS numbers: 71.27.+a, 71.10.Fd, 74.70.Kn, 71.45.Lr 
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I. INTRODUCTION 
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The competition between charge ordered, metallic, and superconducting phases is relevant to a broad range of 
strongly correlated electron materials. For example, in the vanadium bronze /3-Nao.33V20s, superconductivity ap- 
y ■ pears close to a charge ordered phase under an applied external pressure*. The appearance of a pseudogap in 
oxygenated samples of Ndi.ss Ce .i5 Cu 4+a has been suggested to be due to charge ordering*. Quarter-filled layered 
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organic materials such as the BEDT-TTF family of organic molecular crystals 3 with the 6 and (3" molecular stack- 



ing patterns also display a subtle competition of metallic, insulating, charge ordered, and superconducting phases*. 
Superconductivity in organic compounds is usually found in close proximity to ordered insulating phases*. For ex- 
ample, k-(BEDT-TTF)2Cu[N(CN)2]C1 is an antiferromagnetic Mott insulator which becomes superconducting under 
pressure*. Superconducting #-(ET) 2 l3 and /3"-(BEDT-TTF)3Cl2(H 2 0)2 are close to charge ordered phases*. Super- 
conductivity occurs in the quasi-one-dimensional Bechgaard salts, TMTSF2X, when a spin-density wave is suppressed. 
It is then important to understand the connection of the superconducting state to the nearby ordered phases and 
analyze the effect of the fluctuations associated with the ordering transition on the normal metallic phase. 

Several anomalous properties have been observed close to the charge ordering transition in quarter-filled organic 
conductors: (i) Suppression of Drude weight and enhancement of optical spectra at low frequencies at about 500-1000 
cm -1 in metallic 9 3| /3" |j and a-salts |10| at low temperatures, (ii) The temperature dependence of the resistivity 
may be different from Fermi liquid behavior, in particular, the resistivity can increase as the temperature is decreased 
just before becoming superconducting (see the Table in Ref. 6). Previously we have explored, using slave bosons, the 
possibility of superconductivity 6 and the metal-insulator transition 11 in the quarter-filled extended Hubbard model. 
Here, we concentrate on the dynamical properties in the metallic phase close to the charge ordering transition. Wc 
find that due to the scattering of electrons from charge fluctuations with (ir, it) wavevector, dynamical and transport 
properties display behavior different from that expected in a typical metal. For instance, a strong suppression of 
quasiparticle weight as well as enhancement of spectral weight at low but finite frequencies takes place as the charge 
ordering transition is approached from the metallic side. Also we examine the possibility of superconductivity mediated 
by short range charge fluctuations close to the transition using both Lanczos diagonalisation and large- N approaches. 
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We find that superconductivity with d xy symmetry is possible close to the charge ordering transition. We note that 
the present analysis is similar in spirit to those that aim to understand the effect of spin fluctuations on the metallic 
phase and the possibility of superconductivity mediated by thern^ in high-T c compounds, k-(BEDT-TTF)2X [lflj . 
heavy fermionSfi* and ruthenatesi*. 

The paper is organized as follows. In Section [HI we introduce an extended Hubbard model used to describe the 
electronic properties of quarter-filled layered molecular crystals. We also review the path integral formalism written in 
terms of Hubbard operators and the large- TV expansion introduced to compute electronic properties of the model. In 
Section IlIII we show results for the dynamical charge correlation function, spectral density, and optical conductivity 
computed with Lanczos diagonalisation comparing them with large- TV results. In Section IV we discuss our results 
contrasting them with available experimental data on the quarter-filled organics. Section V is devoted to the possibility 
of having superconductivity in the model. 

II. DYNAMICAL PROPERTIES IN THE U-INFINITE LIMIT: LARGE- TV APPROACH 

We consider an extended Hubbard model at one-quarter filling on a square lattice. This has been argued to be the 
simplest model needed to understand the electronic properties of the layered molecular crystals with the 9 and (3" 
molecular arrangements within each layer ^ The Hamiltonian is 

H = —t ( c L c j> + c ] CT c«0 + U^2n^mi 

<ij>,<7 i 

+ v y] mnj - ny^jg (l) 

<ij> in 

where U and V are the on-site and the nearest-neighbors Coulomb repulsion, respectively. c\ a creates an electron of 
spin a at site i. In the limit U 3> V 3> t the ground state is insulating with a checkerboard charge ordered pattern*. 
For U — > oo and V = 0, the system is expected to be metallic as it is quarter-filled. Evaluation of the Drude weight 
by Lanczos techniques suggests a metal-insulator transition takes place at a finite value of V c ss 2.2t for a sufficiently 
large valueii of U = lOt. 

We now introduce the Hubbard operators 

X°° = {l-c\- a c m )c ia , Xf = {X°°)\ xr'=4^>. (2) 

The five Hubbard AVoperators Xf a and Xf° are boson-like and the four Hubbard A-operators Xf and X® a are 
fermion-like. The names fermion-like and boson-like come from the fact that Hubbard operators do not verify the 
usual fermionic and bosonic commutation relations^*. 

In the [/-infinite limit, the Hamiltonian (t — J — V model with J = 0) can be written in terms of Hubbard 
operators as 

H{X)= J2 t ij X?°X9°+ ]T VijXrxf-n^Xr- (3) 

< ij > , u < i j > , tr i,a 

where fi is the chemical potential. The Hubbard operators in this limit satisfy the completeness condition 

x™+j2xr = i, (4) 

a 

which is equivalent to imposing that "double occupancy" at each site is forbidden. 

There are two main difficulties in the calculation of physical quantities using Hamiltonian QJ: the complicated 
commutation rules of the Hubbard operators^* and that there is no small parameter in the model. A popular method 
for handling the former difficulty is to use slave particles. For instance, within the slave boson method^, the original 
fermionic X 0(T operator is decoupled as X 0<7 — b'f a> where b and / are usual boson and fcrmion operators, respectively. 
The second difficulty can be dealt with by using a non-perturbative technique (which we will use in the present paper) 
based on a large- N expansion, where N is the number of electronic degrees of freedom per site and 1/N is assumed 
to be a small parameter. At one-quarter filling (which is the main interest in this paper), we expect the large- TV 
approach to be a good approximation. This has been shown in the overdoped regime of high-T c cupratesi^. 

Hamiltonian has been treated via large- TV in a slave boson representation in Ref. [IT] for V = 0, and in the 
context of quarter-filled layered organic superconductors (V ^ 0) in Ref. |j. Here, we concentrate on the dynamical 
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properties of Hamiltonian (JTJ, using the recently developed large- A expansion^. This method is based on a path 
integral representation of the Hubbard A-operators which is written in terms of Grassmann and usual bosonic variables 
associated with fermi-like and boson-like operators, respectively. In doing this, additional constraints are needed to 
make these field variables behave as Hubbard operators (satisfying their associated algebra), as they should. Although 
this may seem a great complication in the theory, in fact it avoids introducing any decoupling scheme of the original 
Hubbard operators, as in slave boson representations. For completeness we will summarize the framework used in the 
diagrammatic expansion developed in Ref. fl9| . 

Our starting point is the partition function Z written in the Euclidean form 



crO vOa' 
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X90 

x (sdetM AB )?ex ] » (- / ,/r L L [X.X)) . (.".I 
The Euclidean Lagrangian Le(X,X) in J5J is 

L E (X, X)= 1 -J2 X °° + X ^ + H(X) . (6) 



i.a 



It is worth noting at this point that the path integral representation of the partition function [5J looks different to 

that usually found in other solid state problems. The measure of the integral contains additional constraints as well 

i 

as a superdeterminant, (sdetM^s)^ . Also the kinetic term of the Lagrangian (JSJ) is non-polynomial. The determinant 
reads 

(sdetA/ AB )? = (7) 

i 

and is formed by all the constraints of the theory. Note that (sdetAf^s) ? is not proportional to (—A 00 ) 2 , because the 
theory is constrained in a supersymmetric sense where boson and fermion determinants must be treated in different 
ways (see Refpjl for more details about the path integral formalism for Hubbard operators). The constraints appearing 
in the theory are necessary in order to recover the correct algebra of the original Hubbard operators. In Eq. <|10[) we 
show how to treat this determinant through the use of a large- A expansion. 

We now discuss the main steps needed to introduce a large- A expansion of the partition function (JSJ. First, we 
integrate over the boson variables X aa using the second ^-function in |JSJ. We extend the spin index a = ±, to a new 
index p running from 1 to A. In order to get a finite theory in the A — > oo limit, we re-scale the hopping £y to Uj/N 
and Vij to Vij/N. In doing so, note that Uj/N (rather than ty), should be fitted to band structure calculations. The 
completeness condition is enforced by exponentiating: Xf° + J2 P Xf = N/2, with the help of Lagrangian multipliers 
Xi. We write the boson fields in terms of static mean-field values, (ro, Ao) and dynamic fluctuations 

A,™ = Nroil + SR,) 
A, = Xo + SX u (8) 

and, finally, we make the following change of variables 



Jip VNF * 

f ip = -J= A 0p , (9) 
VAr 

where and fi P are Grassmann variables. 

Introducing the above change of variables (Eq. (JSJ and Eq. iJSJ) into Eq. JBJ and, after expanding the denominator 
appearing in ©, we arrive at the following effective Lagrangian: 

iV 

2 



i,p 

N N 

+ E ^ r oftpfjP + E Vi^SRiSRj -fxJ2 ftpfivi 1 ~ SRi + 5Rj 



i,J,P i,3 hP 
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i i,p 

- £ Z\ p (1 - 5R t + 5Rj) Z ip , (10) 

ip 

where Ao has been absorbed in the chemical potential /i — > /i — Ao and all constant and linear terms in the fields have 
been dropped. The path integral representation of (sdetM^B)^ ; written in terms of the iV-component boson ghost 
fields^ Zp, leads to the last term of Lagrangian illUfl . Note that all the complications arising from the Hubbard 
algebra have been translated to an effective theory of fermions interacting with bosons. Indeed, the interaction terms 
appearing in the effective lagrangian (|10|l are generated solely by the Hubbard algebra (apart from the no double 
occupancy constraint) and are not present in the original hamiltonian J2J), which is quadratic in the Hubbard operators. 

In the above expansion we have only retained the first non-trivial terms that couple the fermionic and bosonic 
modes. In order to have a systematic scheme to classify and deal with these interaction terms we introduce a set 
of Feynman rules in powers of l/i\Si*. These will help us to determine, for instance, that the terms retained in the 
effective Lagrangian (|1(J[) correspond to expanding through 0(1/N) in the large- N expansion. The Feynman rules 
needed to carry out this project can be summarized as follows: 

(i) Propagators: We associate with the two component 5X a = (6R , SX) boson field, the bare propagator Dq 

4Vr%(cos(q x ) +cos(q y )) r 



which is represented by a dashed line in Fig. ^ connecting two generic components a and b. q and v n are the 
momentum and the Matsubara frequency of the boson fields, respectively. 
The bare propagator of the A^-component fermion field f p reads 

G(ow(k, LO n ) = —- y- (12) 

iu> n - (£k - A*) 

which is represented by a solid line in Fig^ connecting two generic components p and p' . The electron dispersion 
relation appearing in Eq. (|12f) is the one associated with the original fermions renormalized by the interaction: 
Sk = —2tr (cos(k x ) + cos(fcj,)), with t the hopping between nearest neighbors sites on the square lattice. The 
quantities k and w„ are the momentum and the fermionic Matsubara frequencies of the fermion field, respectively. 
We associate with the TV-component bare ghost field Z p the propagator 

T-^pp' = o~pp f (1^) 
which is represented by a dotted line in Fig. ^connecting two generic components p and p'. 

(ii) Vertices: The expressions of the different three- leg and four- leg vertices are 

Af' = - + </) + /*; l) 5 pp ' (14) 
representing the interaction between two fermions and one boson (see Fig^a)); 

a pp' = _f -f (w + - M -| \ 5 pp'^ (15) 

representing the interaction between two fermions and two bosons (see Fig. |%{ b)); 

r; pl = (-i)(v . o), (i6) 

representing the interaction between two ghosts and one boson (Fig. ^c)); and 



— ( 1) ( Q Q ) <W 



(17) 



representing the interaction between two bosons and two ghosts (Fig. ^d)). Each vertex conserves momentum and 
energy, as it should. 
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(a) (b) (c) (d) 




FIG. 1: Feynman diagrams in the large- N expansion of the Hubbard operator theory. Solid lines represent fermions which are 
related to the electrons. Dashed lines represent bosons which are related to charge fluctuations. Dotted lines represent ghosts 
which are not physical but related to the constraints appearing in the theory which enforce that fermions satisfy the Hubbard 
operator algebra, (a) to (d) show the types of vertex which occur up to order 0(1/N). (a) is the vertex between two fermions 
and one boson, (b) is the vertex between two fermions and two bosons, (c) is the vertex between two ghosts and one boson, 
(d) shows the vertex between two bosons and two ghosts, (e) represents the sum of all one-loop diagrams contributing to the 
irreducible boson self energy which is of 0(1/N). 

In the lowest order of the expansion N = oo, we have the original fermions renormalized by the interaction, 
£k = — 2tr (cos(k x ) + cos(ky)). For a given value of fj,, r$ must be determined self-consistently. For instance, tq is 
equal to 5/2 (where S is the hole doping away from half-filling) from © and the completeness condition Q. 

The path integral JSJ is written in terms of the original X-operators without having to introduce slave particles. 
Eq.© is analogue to the path integral used for the Heisenberg model where, using SU(2) coherent states, the measure 
can be written^! in terms of the spin S. There is, however, an extra price we have to pay if we work with the original 
Hubbard operators. For instance, we need to introduce a new constraint (the second delta function in Eq. (J5J) and 
the determinant (sdetAf^s) of the matrix formed by the constraints appearing in the theory 19 . In spite of these 
"apparent" complications our formulation is very flexible in calculating the physical quantities of interest, as it will 
be shown below. 

In summary, we have developed a diagrammatic technique appropriate for a large- N expansion along the lines of 
the large- N expansion developed in quantum field theory. Hence, from the order of the propagators and vertices, we 
can determine the order of the diagram contribution. 

To conclude this section we make contact with closely related approaches such as slave boson formulations. In 
contrast to slave boson theories: (a) Greens functions are calculated in terms of the original Hubbard operators, (b) 
fermions, fi P , appearing in the theory are proportional to the Fermi- like X-operator X op (see to all orders in the 
1 /N expansion; not only to leading order—, (c) as our path integral is written in terms of X-operators we do not need 
to introduce a priori any decoupling scheme, and (d) ro is the mean value of X 00 which is a real field associated with 
the number of holes (see Eq. JSJ) and not with the number of holons. At leading order (N — * oo or O(l)) and V = 0, 
our formalism is equivalent to slave boson approaches. However, at the next to leading order (0(1/N)), (which 
is necessary to calculate one-electron properties such as the electron self energy S(k, u>) and the electron spectral 
function A(k, ui)), the two formulations do not coincide. The differences between the two formulations are not yet 
completely established. Our theory has the significant advantage that it does not require the introduction of gauge 
fields like in slave boson approaches. Hence, through order 0(1/N) we do not need to take care of gauge fluctuations 
nor Bose condensation (note that Eq. [S]does not mean Bose condensation). This is important because for the doped 
Hubbard model the gauge fluctuations are known to significantly change the physics 23 . Careful numerical work will 
determine the improvements of the present approach with respect to slave boson formulations. 

III. DYNAMICAL PROPERTIES OF THE METALLIC PHASE CLOSE TO THE CHARGE ORDERING 

TRANSITION 

In this section, we analyze using large- N and Lanczos techniques the influence of the charge ordering transition on 
the dynamical properties of the normal metallic phase. 
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A. Charge response 



The dynamical electronic density-density response function can be written in terms of Hubbard operators. We 
define the retarded density-density, Green's function as 



f 

r 

pq 

From ^2 q Xf 1 = N/2 — and J3J) we find, after Fourier transforming 



b *i W = N ^ < r r^r(T)Xf (0) > . (18) 



r 2 

D(q, i/ n ) = -iV(-) D M (q, i*,). (19) 

Here Z?fli?(q, is the (i?, i?) component of the boson propagator. This is the only physical component of the boson 
propagator and encodes the charge fluctuations occurring in the system. Other components of the boson propagator 
such as the (A, R) or (A, A) contain the nonphysical field A which are introduced to enforce the no double-occupancy 
constraint. Unlike in slave boson theories, the (R, R) component used here is associated directly with the charge and 
not with a fictitious bosonic field (holon). 

Through 0(1 /N) the boson propagator consists of the bare boson propagator Dm) (which is of order 0(1/ N)) 
renormalized by a RPA-type series of electronic bubbles. The irreducible boson self-energy components, H a b, are 
obtained (through order 1 /N) from the summation of all the contributions corresponding to the one-loop diagrams 
shown in Fig. Qfe). 

The last two diagrams appearing in Fig. ^e) involving ghost fields are very important. It is possible to show that 
these two diagrams exactly cancel the infinities, due to the frequency dependence of our vertices, of the two first 
diagrams appearing in Fig. H^e). Ghost fields interact only with the boson fields as can be seen from Fig. ^c) and 
Fig. Old). Summarizing, the only role of ghost fields, through order 1/N, is to cancel infinities in the boson self-energy 
n a b arising from the frequency dependence of our vertices l|14|) and (|15fl • 

Using our Feynman rules, we can now write out explicitly each of the components of the boson self-energy H ab 

ttRBiq^n) = -tt [2 n F (e k - /i)(e k+q - e k ) 



N s 4 

k 



/ , ,2 [«F(£k+q - M) - «F(£k - M)] 
(£k+q + £k) 



-lV n + £ k + q — £k 



(20) 



IWq,z/„) = -jj- 7)2_^ (£k+ q + £k) 



N s 2 

6 k 



[rt F (g k+q -fi)- ?if(£ k - M)] / 21 s 

— W n + £ k +q — £ k 



and, 



tt / \ w V [»F(£k+g ~ lA ~ n F (sk - //) , 

IIaa q, v n ) = -— > ; ■ . (22) 

N s -%v n + e k+q - e k 

where N s is the number of sites of the system. 

From Dyson's equation and n a ;, the dressed components of the boson propagator, D a b, can be found: 

(D ab )- 1 = (D (0)ab )- 1 -IL ab (23) 

D ab may contain collective excitations such as zero sound^. 

In order to look at charge ordering instabilities induced by the intersite Coulomb interaction, V, we have calculated 
the static charge susceptibility D(q,u n = 0) for different q vectors on the Brillouin zone (BZ). At one-quarter filling 
(S = 0.5) the corresponding chemical potential is fi = — 0.360t in the limit N — > oo. We find that the static 
susceptibility diverges at the wavevector q c = (tt, tt) for V — V c « 0.65i signalling the instability to a checkerboard 
charge density wave (CDW). The value of V c is slightly smaller than the one previously found using slave bosons^, 
V c w 0.69i. This is because of the decoupling of the electron operators introduced within slave bosons to treat the 
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intersite interaction term, Vriirij which is not needed (due to the use of Hubbard operators) in the present large- N 
approach. For comparison, recent exact diagonalization calculations^ give a critical value for the metal-insulator 
transition driven by V at about V c f» 2t for U — 20t. The large difference in V c between large- N and Lanczos 
diagonalization calculations can be attributed to the strong renormalization of the bare band (given by vq = 5/2) 
which appears in large- TV approaches at 0(1). Introducing higher order terms in the l/7V-expansion may give larger 
values of V c , in closer agreement with Lanczos calculations. 

In Fig. |21 we show the evolution of — Im.D(q c , u) as the system is driven close to the charge ordering instability, 
V < V c , for the wavevector: q c = (7r,7r). The intersite Coulomb repulsion softens the plasmon mode at q c which 
appears for U — > oo and V = and, at the same time, increases its weight. At wavevectors far from q c the plasmon 
mode shows up as a peak located at frequencies of about t which carries small weight and is barely influenced by 
V . Because the mode at (tt, tt) is plasmon-like it can be detected, in principle, with electron energy loss scattering 
(EELS) 24 or inelastic X-ray scattering 2 ^. With EELS one is able to obtain information on the electronic properties 
of the system at a given energy and wavevector, so that, for instance, the dispersion relation of the mode can be 
mapped out. A more appropriate way of detecting the plasmon-like mode is by using inelastic X-ray scattering, which 
provides a direct probe of the dynamical charge correlation function and has been succesfully applied to one- and 
two-dimensional Mott- Hubbard systems 2 ^. 

In order to compare with large- TV we compute, with Lanczos diagonalization, the spectral decomposition of the 
charge correlation function 

C(q,u) =Y,\HN tl \0)\' 2 5(v-(E m -E )) (24) 

m 

where iV q = l/vi£; e jqRi (c^Ci — (cfci)). E m and Eq denote the excited and ground state energies of the system, 
respectively. L is the number of sites in the cluster. Note that C(q, v) can be compared to — Im£)(q c , v) as they have 
equivalent definitions. Of course, attention must be paid to the fact that we are comparing calculations of the charge 
susceptibility on an infinite system with calculations on a L = 16 cluster. Indeed, we find that C((\ c ,v) is in rather 
good agreement with — ImD(q c ,is) (see Fig. |2J), both displaying similar softening and increase in amplitude of the 
plasmon mode at (tt, tt) close to the charge ordering transition. 

The imaginary part of the charge correlation close to the charge ordering wavevector: q — > q c can be fitted to the 
following RPA formSi 

-lmD(d,u)=A V (25) 
v + W q 

where w q = lu + C(q — q c ) 2 , where A and C are constants. ui gives the position of the peak appearing in the charge 
correlation function at (tt, tt) for different V's and goes to zero as V — > V c , measuring the proximity of the system to 
the charge ordering transition (see inset of Fig. We note that the overall behavior of the charge susceptibility is 
analogous to the one of the spin susceptibility in nearly antiferromagnetic metalsSi 2 ^. 



B. Spectral densities 

In order to calculate spectral densities, we first discuss the evaluation of the self-energy corrections to the bare 
fermion propagator (|12fl . which occur at 0(1/ N). Using our Feynman rules there are two diagrams, shown in Fig. 
EU&)i contributing to the self energy to 0(1/N). The analytical expression for these two diagrams reads: 

S PP = E K P 'D ab G p ^A( p + J2K P b D ab (26) 

p\p" ,a.b a, b 

where integration over internal momenta and sum over Matsubara frequencies is assumed. 

Using the spectral representation for the boson fields, D ab , the imaginary part of the self-energy JmS can be 
obtained 

ImS(k,w) = -^E{ilm( J D M (q,w-e k _ q ))(e fc _ 9 + 2/i + w) 2 
s q 

+ lm(D R x(q,oj - e k _q))(£ k _ q + 2fi + u) 

+ Im(£)Ax(q,w - £k- q ))}(ris(w - £k-q) +rii?(-£k-q)) (27) 
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FIG. 2: The softening of the plasmon-like mode at the wavevector q c = (ir, n) as the system is driven closer to the checkerboard 
charge ordering transition. The frequency dependence of the charge correlation function is shown for several different values 
of V/t. The right and left panel show results obtained using Lanczos diagonalization on L — 16 site clusters (U = 20t) and 
large- N theory, respectively. A Lorentzian broadening of r\ — O.li is used in the calculations. Only for wavevectors close to 
or at (ir,n), the softening of the plasmon mode is observed as a consequence of the proximity of the system to a checkerboard 
charge ordering transition. Calculations of dynamical properties using large- N theory at 0(1/N), which couples the electrons 
to the short range charge fluctuations associated with this transition, and Lanczos diagonalization, suggest that this plasmon 
mode is responsible for the 'unconventional' behavior of dynamical properties. The inset compares the position of the plasmon 
peak at (n, n) computed from Lanczos and large- TV approaches. 



(a) 



(b) 

— r 

i 



FIG. 3: (a) Contributions to the electron self-energy, E(k, ui), through 0{1/N), in the Hubbard operator theory. The first 
diagram contains two three-leg vertices as the ones shown in Fig. 0]a) and the second one is formed with one four-leg vertex 
as shown in Fig. 0^b). (b) Contribution to the effective interaction between quasiparticles, V e //, through 0(1/N). This 
interaction is used in the present work to analyze superconducting instabilities of the Fermi liquid induced by the charge 
fluctuations appearing close to a checkerboard charge ordering transition induce by V. 




FIG. 4: Evolution of the real and imaginary parts of the self-energy of an electron at the Fermi surface as the system is driven 
close to the checkerboard charge ordering transition from large-iV theory. The amplitude of the self-energy is enhanced at 
frequencies between t and 3t due to the enhancement of fluctuations associated with (it, tt) short range charge ordering. The 
behavior of the self-energy leads to an enhancement of spectral weight in the spectral density (see Figs. |5J and an incoherent 
band in the DOS (see Fig. |HJ between u> = t and 3t as we approach the charge ordering transition. The intersection of the 
curve of ReE(w) versus u> with u + n — e(k) determines the quasi-particle peaks in the electronic spectral function. 



Note that this self-energy is the one associated with the propagator G(k, ui) of the Hubbard X-operators. In 
contrast, in slave-boson approaches a convolution of the fermion and boson operators is needed in order to recover 
the actual electronic self-energjiSS. 

Fig. QJshows the behavior of ImE(k, oj) with increasing V for a wavevector on the Fermi surface: k = (1.204, 1.204) 
(we have used r] = O.lt in the analytical continuation). As we approach V = V c , both the imaginary and real parts of 
the self-energy (which from Eq. I|27l involves a sum over the full BZ) are enhanced in the positive range of frequencies 
t — 3t due to the scattering of the electrons off the checkerboard charge fluctuations. Performing a Kramers-Kronig 
transformation on ImE, we can obtain the real part of the self-energy, ReE, which is also plotted in Fig. 0] 

From ImE and ReE, we can compute the electron spectral function A(\c,u>) = — -ImG(k, ui) as 

. . 1 ImE(k, ui) . 

( ,W > = ~^(cj + /i-£ k -ReE(k,tj)) 2 +lmE(k,w) 2 ( ' 

In Fig. Elwe show the spectral function obtained from large- TV theory, for an electron at k = (0, 0), (1.204, 1.204), (tt, 7t), 
for different values of V — > V c . The spectral density of an electron at the Fermi surface displays a quasiparticle peak 
characteristic of a Fermi liquid at ui = ji. The rest of spectral weight that is left is incoherent. 
The quasiparticle weight, Z^, evaluated at the Fermi surface is defined as 

Zk = (1 _^|l^ rVo . (29) 

In the inset of Fig. [5] we observe how a gradual suppression of occurs as the charge ordering transition is 
approached. This can be compared to the suppression of the Drude weight found in Lanczos calculations^, which is 
also evident in the spectral function plotted in Figs. El and Spectral weight is transferred from the quasiparticle 
peak to the range of energies between t — 3t, as V tends to V c due to the scattering of the electrons off the charge 
fluctuations associated with short range checkerboard charge ordering. The modes close to (tt, tt) give the strongest 
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contribution to the scattering. The apparent peak around uj — — 2t should not be interpreted as a quasiparticle peak 
but as the lower Hubbard band^ associated with the on-site Coulomb repulsion U. 

The behavior of the spectral density shown in Fig. can be now understood from the evolution of the real part of 
the self energy shown in Fig. 0] The scattering of electrons from the strong charge fluctuations at (jr, it) wavevectors 
involves large frequencies. This leads to an enhancement of the real part of the self-energy at large frequencies which, 
in turn, produces an increase of spectral weight at large and intermediate energies. This behavior is analogous to the 
one found in metals in the presence of short range spin fluctuations*^. Unlike in the case of nearly antiferromagnetic 
metals, no new poles induced by the interaction arise. This is because at quarter-filling no two points of the Fermi 
surface are connected by the scattering wavevector q c = (it, tt), and therefore the effect of the fluctuations near q c on 
the electrons is weaker than spin fluctuations in systems close to half-filling. 

In order to test the validity of the large- N approach we have also computed the spectral densities from Lanczos 
diagonalization of finite clusters^ 

AW(k,w) = ^|(™,A e + l| c + CT |0,Ay| 2 <5( w - (E m (N e + l)-E (N e ))) (30) 

for adding and electron to the system with N e electrons and 

i4<->(k,w) = £ |<ro,JV e - l|c kCT |0, N e )\ 2 6(tu + (E m (N e - 1) - E Q (N e ))) (31) 

m 

for removing an electron from the N e electron system. E m and E$ denote the excited and ground state energies of 
the system and cjf =1/ VL J^j e lkRj c] CT . 

In Figs. E]and[3we plot the evolution of the spectral densities calculated with Lanczos techniques for wavevectors 
at k — (0, 0), (n/2, 0), (n/2, tt/2) and {tt,tt) for different values of V/t. At k = (0,0) two sharp peaks are clearly 
distinguished already for V = 0. One of them is the quasiparticle peak and we associate the lower one with the lower 
Hubbard band due to the presence of the U. For the nearest wavevectors to the Fermi energy: k = (jr/2,0) and 
(7r/2, 7r/2), we find an enhancement of incoherent spectral weight at finite frequencies as the charge ordering transition 
is approached. 

Finally, the total density of states (DOS) can be computed from 

N(u) = - J2( Ai ~ ] ( k < ^ + Ai+) ( k - ( 32 ) 

k 

In Fig. IHJwe compare the evolution of the DOS, N(ui), for increasing V/t, calculated with both Lanczos at U — ► oo 
and large- N. 

From Lanczos calculations we observe (left panel of Fig. |8J for V = a band at about — 3t, a quasiparticle band 
at situated at u> = fj, and a band running from t to 5t. As V/t is increased the weight of the quasiparticle peak is 
reduced and weight between 2t and 5t is gradually enhanced. Also a suppression of spectral weight at low frequencies 
occurs as a precusor effect before the charge ordering transition takes place. This general behavior is in qualitative 
agreement with large-A^ calculations. Indeed, an incoherent band at negative frequencies of about — 2t, associated 
with the lower Hubbard band, a suppression of states close to the Fermi energy and an overall enhancement of spectra 
between t and 3t occurs (see right panel of Fig. [BJ. However, we note that the pseudogap appearing within large- N is 
less pronounced than in Lanczos calculations. This can be attributed to finite size effects appearing in small cluster 
Lanczos calculations. 



C. Optical conductivity 

It is interesting to analyze the behavior of the optical conductivity as the system is driven through the charge 
ordering transition. Using Lanczos diagonalisation we have computed 

f ,\ nxt \ , 7re 2 \ - \(n\j x \0)\ 2 x/ . , . 

a(u) = DS{u) + — > — — d{u - E n + E ), (33) 

where j x is the current in the x-direction, Eq the ground state energy and E n the excited energies of the system, e is 
the electron charge and the Drude weight is denoted by D. 
The following sum rule^ is satisfied by <t(lu), 

r°° -ne 2 
J cr{uj)doj = -— < 0\T\0 > . (34) 
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FIG. 5: Evolution of quasiparticle spectral density of states computed from large- N theory of an electron at the wavevectors 
k=(0,0), (fc_F, &f), and (tt, n). Close to the charge ordering transition spectral weight is transferred from the quasiparticle peak 
to low and intermediate frequencies. The quasiparticle weight at the Fermi energy, Zk, is rapidly suppressed (see inset) as the 
charge ordering transition is approached V — > V c . The results presented here can be compared with A(\f.,u)) computed from 
Lanczos shown in Fig. |S|and[7| 
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FIG. 6: Evolution of quasiparticle spectral density of states, A(k,u>), at k = (0,0) and (0,7r/2), computed from Lanczos 
diagonalisation on a L = 16 cluster for an extended Hubbard model at quarter-filling. The on-site Coulomb repulsion is taken 
to be U = 20t and a broadening of the delta peaks, 77 = O.lt is used. As the system approaches the metal-insulator transition, 
an enhancement of spectral weight at finite frequencies and a suppression of the weight at the Fermi energy takes place. At 
k = (0, 0), the two sharp peaks are associated with the lower Hubbard band and the quasiparticle peak. An overall qualitative 
agreement with the results from large- N theory is found (see Fig. |SJ. 



where T is the kinetic energy operator, which is the first term in the Hamiltonian (1). 

The optical conductivity is plotted in Fig. ^\ for increasing values of the ratio V/t and fixed U = 20t. At V = we 
find a Drude peak and a broad mid-infrared band centered at about 3t. As V/t is increased the mid-infrared band is 
enhanced and a well-defined feature builds up at the lower edge of the mid-infrared band, at frequencies of about 2t. 
Also an incoherent band present at larger energies of the order of U (not shown for clarity) is gradually suppressed and 
its associated weight transferred to the mid-infrared band as V is increased. From the behavior of spectral densities 
and DOS shown in Figs. [5J we attribute the enhancement of optical weight observed in the mid-infrared range 
to an increase in the incoherent excitations carried by each quasiparticle as a result of charge fluctuations associated 
with short range checkerboard charge ordering. From the behavior of the spectral densities, A(k, u>) shown in Fig. 
El and [7| and assuming that a lowest order diagram (neglecting vertex corrections) is enough to compute the optical 
conductivity we would attribute the low energy feature to transitions between the incoherent band carried by each 
quasiparticle and the quasiparticle peak situated at the Fermi energy. This interpretation is plausible if one notes 
that the low energy feature observed in Fig. Amoves together with the broad band as V/t is increased. Similar results 
would be obtained from large- N theory evaluating the bubble Feynman diagram for the optical conductivity, as the 
spectral densities obtained are similar to the ones obtained from Lanczos diagonalization. 



IV. CONNECTION TO EXPERIMENTAL RESULTS 



Recent experiments with Raman scattering^ and optical conductivity measurements* on the insulating salt 6- 
(BDT-TTP)2Cu(NCS)2 find that the checkerboard charge ordered state discussed in this paper is indeed the ground 
state. This gives experimental support to the model discussed here. A discussion of other possible orderings within 
more complicated models can be found in the work by Seo^i and Clay, Mazumdar, and Campbell^. 

We review now the experiments on resistivity measurements on several quarter-filled organics, and make contact 
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FIG. 7: Evolution of quasiparticle spectral density of states, A(\l,u), at k = (n/2,n/2), and (tt,tt), computed from Lanczos 
diagonalization on a L = 16 site cluster for an extended Hubbard model at quarter-filling. The on-site Coulomb repulsion is 
taken to be U = 20t and a Lorentzian broadening of the delta peaks, rj = O.lt is used. 




FIG. 8: Evolution of the total density of states (DOS) as the charge-ordering transition is approached from the metallic phase. 
The left and right panels show results from exact diagonalisation on a 16 site lattice with U — > oo and large-iV approaches, 
respectively. (The critical value of V is V c ~ 2i and 0.65i, respectively). A Lorentzian broadening of rj = O.lt has been 
introduced in the exact diagonalisation calculations, to aid comparison with the large N results. As the intersite Coulomb 
repulsion V is increased, the density of states close to the Fermi energy is gradually suppressed indicating the proximity to the 
charge ordering transition. At the same time spectral weight is enhanced for frequencies in the range t to St in the large- TV. 
The peak at -2t is an incoherent band associated with the lower Hubbard band. An overall qualitative agreement between 
Lanczos and large- N calculations is found. 
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FIG. 9: Evolution of the optical conductivity computed from Lanczos diagonalisation as the system is driven through the 
metal-insulator transition. The calculation is performed on a L — 16 site cluster, U = 20t and different V , with a Lorentzian 
broadening of r\ = OAt. Enhancement of optical weight at low frequencies is found as V is increased. The broad band situated 
at about 3V (for large V/t) is due to incoherent transitions between different sites induced by the intersite Coulomb repulsion. 
We interpret the low energy feature appearing at about 2t as a consequence of transitions between the incoherent band and 
the quasiparticle peak found in the spectral densities A(k,ui) for wavevectors on the Fermi surface. 



with the predictions of the large-TV approach presented. 

From the imaginary part of the self-energy (Eq. (J23) we can obtain the behavior of the inverse of the lifetime 
of the quasiparticles, l/r(T) = — 2Im£(fcF,0) with temperature as shown in Fig. ^| From this plot we obtain a 
temperature scale, T* sa 0.22i, at which 1/t(T), changes from T 2 to T behavior. The temperature scale defined 
by T* decreases only slightly as we approach the charge ordering transition remaining finite as V — > V c (through 
0(l/iV). This is in contrast to dynamical mean-field approaches where a similar low temperature scale is suppressed 
as the Mott-Hubbard metal insulator (driven by U instead of V) is approached 37 . Hence, Fermi liquid behavior is 
found below this temperature scale even close to the charge ordering transition occurring at V ~ V c . Presumably, 
higher order corrections in the 1/N expansion may suppress the region where the system behaves as a Fermi liquid 
as V — » V c . Future work should focus in understanding this issue better. 

We have also computed the temperature dependence of the effective mass defined as: m*/m = l/Zk(T), evaluated 
at the Fermi surface and is shown in Fig. 1111 Large- N theory predicts an increase of m* jra as the temperature is 
raised for V — > V c . This means that the system becomes more incoherent as the temperature is increased. Interestingly 
this behavior is also found in the Hubbard model in the limit of infinite dimensions close to the Mott metal-insulator 
transition 3 ®. However in that case the system is close to a metal-insulator transition which takes place between two 
non-ordered phases, in contrast to the charge ordering transition discussed here. At V = the effective mass is 
temperature independent as one would expect from a weakly interacting system. At the lowest temperatures we 
obtain enhanced effective masses in the range 1.3 to 2, for V/t, varying from to 0.63. 

In Fig. [Hwe show optical conductivity data of 0-(BEDT-TTF) 2 CsCo(SCN) 4 along the a- direction^ at T = 15 
K. This salt is metallic down to temperatures of about 10 K. At this temperature a charge ordering transition to an 
insulating phase takes place. The observed optical conductivity displays a band situated at 1.2 eV, a mid-infrared 
band appearing at frequencies of about 0.25 eV, and a feature appearing at low frequencies of about 0.13 eV. In the 
same figure we present a comparison of our exact diagonalization calculations of the optical conductivity performed 
on a 16 site cluster, U/t = 20 and V/t = 1.2 where we have set the hopping energy scale to t = 0.0616^, so that 
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FIG. 10: Temperature dependence of the scattering rate l/r(T) = — 2ImS(kF,^ = 0). This scattering, which is due to charge 
fluctuations, increases as the charge-ordering transition is approached. As the temperature increases above about T* w 0.22t, 
l/r(T) changes from a T 2 dependence to a linear behavior in T. This temperature scale depends only slightly with V, so that 
large- N theory (through 0(1/ N)) predicts Fermi liquid behavior close to the charge ordering transition at temperatures below 
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FIG. 11: Temperature dependence of the effective mass of an electron on the Fermi surface, m*/m = l/Zk(T), as obtained 
from large- AT theory through 0(1/ N). As the system is driven closer to the charge ordering transition a stronger increase of 
the effective mass with T is found. 
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FIG. 12: Comparison of the optical conductivity computed from Lanczos diagonalisation on L = 16 site clusters with ex- 
perimental results— for the metallic salt #-(BEDT-TTF)2CsCo(SCN)4. For the exact diagonalization results we have chosen: 
U — 20t and V — 1.2t. In order to fit the data we chose t = 0.061 eV, which can be compared to values from Hiickel 
band structure calculations^. The lattice parameters for <9-(BEDT-TTF) 2 CsCo(SCN) 4 are a = 9.804A, c = 4.873A and 
Vceii = 4Vmoi = 2074A 3 , where V ce ii and V mo i are the volumes per unit cell and per molecule, respectively. The broad band 
at about 0.25 eV and the low energy feature at 0.13 eV can be explained from short range charge ordering induced by the 
intersite Coulomb repulsion V . This behavior is characteristic of several quarter-filled layered metallic salts which undergo a 
metal-insulator transition at sufficiently low temperatures. 



we associate the mid-infrared band observed experimentally with the one from exact diagonalization calculations. In 
this way, we recover the main features appearing in the experimental data, including the incoherent high-frequency 
feature and the feature appearing at low frequencies. This behavior is commonly observed in metallic 0-salts close 
to the metal-insulator transition 11 and from Fig. 1121 we notice that the low energy feature can be misinterpreted as 
being part of the Drude peak. Caution is in order when comparing our results with experimental data as shown in 
Fig. ^Jbecause some features like the dip appearing at about 0.17 eV have been interpreted in #-(ET) 2 RbZn(SCN)4 
(where a structural transition takes place with lowering temperature) as being caused by the coupling to vibronic 
modes of the ET molecules— More experimental and theoretical work is needed to understand this issue better. 



V. SUPERCONDUCTIVITY 



In the present section we discuss the possibility of having superconductivity close to the charge ordering transition 
induced by the short-range charge fluctuations which appear in the metallic phase. Here we extend the discussion 
presented in Ref.Q and provide full details of the calculations. We also consider the binding energy of holes using a 
Lanczos calculation. 



A. Large-N: pairing symmetry 

Within the large-TV approach, superconductivity is possible at 0(\/N). As we have already seen at O(l), our 
approach describes quasiparticle excitations with renormalized masses. Interaction between these quasiparticles can 
appear at the next-to-leading order of 0(1/N). The effective interaction between electrons are those represented 
diagrammatically in Fig[3Jb); only the three- leg vertex shown in Fig^a) contributes to the effective interaction 
through order 0(1/N). 



17 



6.0 



4.0 - 



cr 

> 



2.0 



0.0 - 



-2.0 - 



-4.0 



V/t=|0 


V/t=0.6 


55 V/t=0.4 






\ / 






\ / V/t=0.2 \ 






v/t=o 











-1.5 



-1.0 



-0.5 



0.0 



0.5 



1.0 



1.5 



FIG. 13: Behavior of the effective potential between quasiparticles, V e f f(qx, q y ), as a function of momenta for increasing V/t 
values along the q x = q y = q direction. As V/t — > (V/t) c , the effective potential becomes negative at the (n, n) points, becoming 
singular at the transition to the checkerboard charge ordered insulator. It is the momentum dependence of the potential shown 
here which leads to the d^ symmetry of the Cooper pairs. This calculation was done using the large-iV approach through 
0(l/N). 



Using the Feynman rules introduced above (see Fig. Efb)), the interaction between the quasiparticles 0(k, k') reads 
6(k- k', uj n - uj n >) = -[n 2 D RR {k - k',w„ - u n <) + 2fj,D RX (k - k',w„ - uj„>) 

+ DAApC-k'.Wn-Wn')] (35) 

where /i is the chemical potential and D a b are the components of the boson propagator which are obtained from 
Dyson's equation 1)23(1. 

In Fig. ^|we plot the dependence of V e //(q = k — k') = 0(q, u> — ► 0)/(<5/2), that is the static limit of the effective 
interaction mediating the possible pairing between the quasiparticles. This clearly shows the development of the 
singularity due to checkerboard charge ordering at the (tt, if) wavevectors. We note that the effective interaction, 0, 
obtained from Ea. H35|) . which is valid on the whole Brillouin zone, coincides only with the one obtained using slave 
boson approaches, when it is evaluated at the Fermi surfaced. This is not true for wavevectors outside the Fermi 
surface. 

In weak coupling, we use this effective potential to compute the effective couplings in the different pairing channels 
or irreducible representations of the order parameter, i (i = (d x 2_ y 2 , d xy , s)). In this way we project out the interaction 
with a certain symmetry. The critical temperature, T c , can then be estimated from: T C i — X.IZujq exp(— 1/|A,|), where 
u>o is a suitable cutoff frequency and Xi are the effective couplings with different symmetries. These are defined a o 17 i 41 : 



A = i jWK|) j(dkyK|)g,(k')y;//(k' - k) gi (k) 

1 (2tt)2 f(dk/\v k \Mkr [ ' 

where the functions: $i(k), encode the different pairing symmetries, and i>k are the quasiparticle velocities at the 
Fermi surface. The integrations are restricted to the Fermi surface. Xi measures the strength of the interaction 
between electrons at the Fermi surface in a given symmetry channel i. If Xi > 0, electrons are repelled. Hence, 
superconductivity is only possible when Xi < 0. In Fig. 1141 we plot the dependence of the effective couplings in the 
possible symmetry channels with V/t. 

We observe that near the charge ordering instability but still in the metallic phase: V < V c , the coupling in the d xy 
channel, Ad x , becomes attractive whereas other couplings become more repulsive. However, we note that the couplings 
are rather small. This implies that critical temperatures are expected to be small. Similar conclusions have been 
reached with large- N treatments of the [/-infinite Hubbard model at V = close to half-filling^. Exact diagonalisation 
results also lead to similar conclusions^ 1 . Although the eigenvalues are small, our results are nontrivial and show a 
tendency of the model to Cooper pairing in the d xy channel. Intuitively one would think that superconductivity is 
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FIG. 15: Schematic plot showing the Fourier transform of the effective potential, Vef f(lx, q y ) for V/t ~ (V/t) c to real space. 
V e ff(x,y) is understood as follows. A quasiparticle is placed at the origin. For instance, if another quasiparticle is placed 
also at the origin there is a large repulsion between them due to the large on-site Coulomb repulsion. This is shown by the 
large positive vertical bar at the origin. At neighboring sites (along the x, and y-directions) the effective potential between 
quasiparticles is also positive, i.e., repulsive. However, at the next-nearest-neighbor sites (along the diagonals of the lattice), 
the potential becomes attractive. This leads to d xy pairing of the quasiparticles. 



less and less favorable when increasing V, due to the repulsion between electrons in neighboring sites. Contrary to 
this intuition we find that short range charge fluctuations can mediate pairing close to V c . 

In Fig. 1131 we observe that as we increase V , the effective interaction becomes more repulsive at small momentum 
transfer. On the other hand, they become more attractive for momenta transfer close to (tt, tt). That d xy symmetry is 
favored can be more clearly understood from Fig. ^] which shows a schematic plot of the Fourier transform of V e f / (q) 
(see Fig. IT5jl . One sees that the potential is negative for an electron placed at the nearest neighbor diagonal sites 
of the lattice while it is positive along the x and y-directions. This is in contrast to the effective potential resulting 
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from spin-fluctuations which show the opposite behavior. This can be understood from previous calculations on a 
3-D extended Hubbard model close to half-filling within RPA performed by Scalapino, Loh and Hirschii, which found 
that the effective potential for charge fluctuations has a negative divergence at (tt, tt) as the transition is approached 
whereas for spin fluctuations it is positively diverging—. Due to the fact that the Fermi surface at one-quarter filling 
is small (no two points in the Fermi surface are connected by the (tt, tt) wavevector), the interaction is less effective 
in inducing pairing as compared to spin fluctuations in nearly antiferromagnetic metals close to half-filling. 

The T c values shown in the phase diagram in Ref. |6( are larger than the values that would be obtained for T C( j XH 
from the BCS equation. In Ref. 6] T c (for each V) was taken to be the temperature below which the coupling Xd x 
becomes negative. Such a calculation is indicative for superconductivity. However, the appropriate way to obtain T c 
is by solving the associated Eliashberg gap equation. 

In conclusion, in the present study we find tendencies to superconductivity in the d xy channel mediated by short 
range charge fluctuations which appear in the metallic phase close to the charge ordering instability. 

B. Lanczos diagonalization: binding energies 

We have computed the binding energy of two holes for different values of V/t and U — 20t on different clusters. 
The binding energy of two holes for L — 16 is defined asM^^ 

E B (2 holes) = (E(6) ~ E{8)) - 2(E(7) - E(8)), (37) 

where E{N e ) is the energy of the system with N e electrons. 

In Fig. 1161 we plot the binding energy for different values of V. Initially, as we increase V, the binding energy 
becomes more positive. This corresponds to the weak coupling regime where one naively expects that V keeps the 
quasiparticles farther apart. Further increasing of V closer to the metal-insulator transition but still in the metallic 
phase leads to a negative binding energy of two holes. From finite size scaling of the binding energy of clusters up 
to L — 20 sites, we find that this happens at about V ~ l.Qt. This finite size scaling is shown in the inset of Fig. 
1161 for values of V close to the metal-insulator transition. Remarkably, the binding energy of two holes changes only 
slightly in the range of values t < V < 2t when going from L = 16 to L = 20. However, in the region V > V c = 2t, 
the results change significantly as we increase the size of the cluster from L = 16 to 20, and the binding energy would 
eventually extrapolate to a positive value in the thermodynamic limit. It is interesting that this region corresponds to 
the insulating phase found earlier— from Lanczos calculations of the Drude weight. An interpretation of our results 
can be made based on previous works^M which studied the binding energy of an extended Hubbard model of the 
high-T c superconductors close to half-filling as a function of V. As V is increased charge fluctuations associated with 
checkerboard charge ordering increase and the quasiparticles existing at small V gradually dress up with a cloud 
of checkerboard charge excitations. This leaves signatures in the one-electron dynamical properties as explained in 
previous sections. Further increasing of V leads to pairing between the quasiparticles mediated by the strong charge 
fluctuations^. Increasing V even further drives the system into the insulating phase. 

Summarizing, a definitive conclusion about superconductivity cannot be made from our results. However, it is 
remarkable that both large- N and Lanczos diagonalization calculations show a similar tendency to pairing of quasi- 
particles in the metallic phase close to the charge ordering transition. Large- N theory singles out the d xy symmetry 
as the preferred pairing channel of the quasiparticles. This symmetry is consistent with the checkerboard charge order 
present in the region where the Lanczos pairing energy becomes negative. 

VI. CONCLUSIONS 

In summary, using a combination of large- N and Lanczos techniques we have explored the dynamical properties of 
the extended Hubbard model at quarter-filling. This is motivated by its relevance to a large class of superconducting 
layered organic molecular crystals. The correlation functions computed from large- N theory and Lanczos techniques 
are found to be in good agreement. Indeed, close to the charge ordering transition driven by the intersite Coulomb 
repulsion, V, several features are found: (i) The quasiparticle weight, Zk, is rapidly suppressed near the charge 
ordering transition, (ii) Spectral density is enhanced at frequencies ranging from t to it, which is also reflected in 
the optical conductivity, (iii) From the computation of the electron scattering rate we find Fermi liquid behavior up 
to T « T* , where T* does not depend strongly on V. For T > T* the scattering rate behaves linearly with T. (iv) 
From large- N calculations we find that superconductivity with d xy symmetry is favored close to the charge ordering 
transition. Exact diagonalization calculations of the binding energy of two holes are consistent with this possibility. 

Given our prediction of unconventional superconductivity in the 9 and (3" molecular crystals it is desirable that 
more measurements be made to test for this. The only evidence so far comes from a measurement of the temperature 
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FIG. 16: Binding energy of two holes for different values of V and U = 20t from exact diagonalisation calculations on L — 16 
and L — 20 clusters. Close to the charge ordering transition, in the range, V > l.&t, binding of two holes becomes favorable. 
V c denotes the value of V at which the metal-insulator transition is estimated to take place from Lanczos calculations of the 
Drude weight. Note that the value of V at which E_b(2 holes) becomes negative is robust against increasing the cluster size 
from L = 16 to L = 20. These results are consistent with large- TV calculations supporting the possibility of pair formation close 
to the charge ordering transition. 

dependence of the London penetration depth of /3"-(BEDT-TTF)2SF 5 CH2CF 2 S03. It was found to go like T 3 at low 
temperatures— This is inconsistent with an s-wave state, but also deviates significantly from the linear temperature 
dependence expected for a d-wave state. On the other hand, the temperature dependence of the heat capacity is 
consistent with s-wave pairing— Electronic Raman scattering could be used to investigate the symmetry of the 
superconducting order parameter. For d xy symmetry, Raman scattering in the superconducting state should show, at 
low frequencies, either to, ui 3 , or lo behavior for Bi g , B^g and, A\ g symmetries, respectively^. 

An important issue to be resolved concerns the role of spin fluctuations in the quarter-filled materials. To first 
order in 1/N, the large- N approach used here does not take spin fluctuations into account^ Measurements of the 
nuclear magnetic resonance relaxation rate and Knight shift should be done in the metallic phase for the relevant 
superconductors. If the spin fluctuations are not important there should be no enhancement of the Korringa ratio. 
This is in contrast to the large enhancements seen in k-(BEDT-TTF)2X superconductors which are close to an 
antiferromagnetic Mott insulator—. 

One way to theoretically investigate the role of the antiferromagnetic spin fluctuations that may be present near 
the charge-ordering transition is as follows. Well into the insulating charge ordered phase (i.e., for V ^> t) it is known 
that there is an antiferromagnetic exchange interaction J' — At 4 /9V 3 that acts along the diagonals of the square 
latticed. Some remnant of this effect will still be present when there is short-range charge order. This could be 
modelled by considering a t — J' — V model where the J' acts only along the diagonals. This model could be studied 
by the same large- N method used previously to study a large family of t — J — V models 51 . There it was found that 
the superexchange, acting along the vertical and horizontal lattice directions, produced d x 2_ y 2 superconductivity. 
Based on that work we anticipate that the effect of the superexchange, which now acts in directions rotated by 45 
degrees, will be to produce d xy superconductivity. Hence, it is possible that charge and spin fluctuations work together 
co-operatively to produce d xy pairing. 
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